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Comment on "Spectral Signatures of the Fulde- 
Ferrell-Larkin-Ovchinnikov Order Parameter in 
One-Dimensional Optical lattices" 

In a recent letter Reza Bakhtiari et al. [1] studied an 
imbalanced two-component atomic Fermi gas in a one- 
dimensional (ID) optical lattice with a trapping poten- 
tial, within the Bogoliubov-de Gennes (BdG) approxima- 
tion. They showed that the prominent oscillations of the 
pairing gap (within the BdG approximation), character- 
istic of a Fulde-FerrcU-Larkin-Ovchinnikov state (FFLO), 
could be detected in the rf spectra and in the momentum- 
resolved photoemission spectra of the gas. In this com- 
ment we show that the BdG approximation not only pro- 
duces inaccurate results for the examples presented in [H , 
but that they are qualitatively incorrect making the anal- 
ysis of the rf spectra unreliable and shedding doubts on 
the applicability of rf spectroscopy to detect the FFLO 
state in ID optical lattices. 

The Letter [l[ is devoted to ID optical lattices with- 
out any reference to higher dimensions where the BdG 
approximation could be valid. The use of the BdG ap- 
proximation is justified by saying that it provides quali- 
tative information on the system and allows one to cal- 
culate the rf spectrum. In spite of the known failure of 
the BdG approximation in ID systems, it is still used 0. 
In the absence of a trap the BCS treatment gives large 
errors in the order parameter as compared to the exact 
solution (see Ref. 19 of the letter). With the inclusion of 
the trap we have to resort to the Density Matrix Renor- 
malization Group (DMRG) to produce numerically "ex- 
act" solutions for ID lattice problems of moderate sizes 
{'-^ 100 sites). Indeed, the DMRG has already been used 
to study the FFLO phase in a ID trapped lattice gas (see 
Ref. 11 and 16 of the letter). Therefore the DMRG con- 
stitutes an ideal benchmark to test the accuracy of the 
BdG approximation in ID lattice problems. 

In Fig. 1 we show the densities and the absolute 
value of the gap for a system with polarization P = 0.23 
(iV| = 40, Ni = 25) in the upper panels and for P = 0.70 
(A^l = 40, Ni = 7) in the lower panels, in a lattice 
of L = 150 sites. Panels A and D correspond to the 
BdG approximation neglecting the Hartree term like in 
Fig. 1 of [H. In panels B and E include the Hartree 
term into the BdG approximation. Panels C and F 
are DMRG results. The gap parameter is defined as 
I Ail = U ^J{< rii-^riii > — < n^i >< >), where the 
expectation values are taken in the ground state wave- 
functions of the three approximations. Within the BdG 
approximations this definition coincides with local pair- 
ing gap. The Hartree term, neglected in [l[, does modify 
the results towards the "exact" DMRG results. How- 
ever, both BdG approximations are still quite far from 
the DMRG results, specially for lower polarizations. The 
most dramatic differences can be seen in the pairing gap 
with very large amplitude oscillations which are com- 
pletely softened in the DMRG results. It might be argued 



FIG. 1: (Color online) Density of majority and minority 
atoms m and absolute value of the local gap parameter |Ai| 
as a function of the site index i. Upper panels show results 
for P = 0.23 and lower panels for P = 0.70. From left to right 
we show results of BdG without the Hartree term, BdG with 
the Hartree term and DMRG. 


that the dimensions of the systems considered are small, 
however, doubling the size of the systems the authors 
found similar results as we did using DMRG. 

The subsequent analysis of the rf spectra is based on 
the BdG results of panels A and D. Taking into ac- 
count that the gap oscillation are reduced by one order 
of magnitude in the DMRG results, it is doubtful that 
rf spectroscopy could provide information about the spa- 
tial structure of the pairing gap. Whether rf spectroscopy 
could signal the FFLO phase in ID lattice systems is still 
an open question which might be confirmed by rf exper- 
iments, or addressed numerically by means of a DMRG 
study. Note that within the Correction Vector Approach 
the calculation of the excitation spectrum is not needed. 
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